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Abstract 

We present the analytical calculation of the four-loop QCD /3-function within the 
minimal subtraction scheme. 



1 



The renormalization group /3-function in Quantum Chromodynamics (QCD) has a history 
of more than 20 years. The calculation of the one-loop /3-function in QCD has lead to the 
discovery of asymptotic freedom in this model and to the establishment of QCD as the theory 
of strong interactions |l]]. 

The two-loop QCD /3-function was calculated in [Q]. Calculations of the three-loop QCD 
/3-function were done in || [| within the minimal subtraction (MS) scheme ||. The MS- 
scheme belongs to the class of massless schemes where the /3-function does not depend on 
masses of the theory and the first two coefficients of the /3-function are scheme- independent. 

In this article we present the analytical four-loop result for the QCD /3-function. Through- 
out the calculations we use dimensional regularization ||] and the MS-scheme. The dimension 
of space-time is defined as D = 4 — 2e, where e is the regularization parameter fixing the 
deviation of the space-time dimension from its physical value 4. 

The Lagrangian for a massless non-abelian Yang-Mills theory with fermions is 



G- u = d^-d^ + gr^Al 

= ^-igA«[T% (1) 

where the gauge-fixing and gauge-compensating parts of the Lagrangian in the covariant 
gauge are 

L gc = d^ a (d^ a -gr bc JA^) (2) 

The fermion fields (the quark fields in QCD) transform as the fundamental representation 
of a compact semi-simple Lie group, 9=1, ...,n/ is the flavour index. The Yang-Mills fields 
(gluons in QCD) A^ transform as the adjoint representation of this group. u) a are the ghost 
fields, and £ is the gauge parameter of the covariant gauge. 

T a are the generators of the fundamental representation and f abc are the structure con- 
stants of the Lie algebra, 

rj-iarpb rjhrjia ^ JdbCrpC 

In the case of QCD we have the Lie group SU(3) but we will perform the calculation for an 
arbitrary compact semi-simple Lie group G. Since the /3-function does not depend on masses 
in the MS-scheme, we will consider the massless theory. 
The definition of the 4-dimensional /3-function is: 

= -(3 a 2 s - Pta 3 s - fcai - 3 a 5 s + 0(a 6 s ) (4) 

in which a s = a s /4-K = <? 2 /167r 2 , g = g(fJ?) is the renormalized strong coupling constant of 
the standard QCD Lagrangian of eq.(Q). \i is the 't Hooft unit of mass, the renormalization 
point in the MS-scheme. 
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To calculate the /3-function we need to calculate the renormalization constant Z as of the 
coupling constant 

as = Z as a s 

where as is the bare (unrenormalized) charge. We obtain this renormalization constant in the 
4-loop order by calculating the following three renormalization constants of the Lagrangian: 
Zhhg for the ghost-ghost-gluon vertex, Zh for the inverted ghost propagator and Z g for the 
inverted gluon propagator. Then Z aa = Z\ hg jZ\jZg. This is from a calculational point 
of view one of the simplest (and most straightforward) ways to obtain this renormalization 
constant at higher orders. However, several other choices are possible such as calculating the 
renormalization factors of the quark propagator, the gluon propagator and the quark-gluon 
vertex. One could also use the background field method f7j, which reduces the calculation of 
the /3-function to the calculation of the gluon propagator only. However, we should note that 
the Feynman rules are more complicated in that case, and this complication would lead at 
the 4-loop level to a complexity of the calculations that is comparable to our more standard 
approach. 

The expression of the /3-function via Z ag is given by the following chain of equations 
d(a B ^ £ ) _ ? 2e i dZ as da s 2e da s 2e 

^ i = ~ fl/^'" 8 = - £a s ~ a sTr~ UsZa}) = ~ £a s + f3{a s ) (5) 

where one uses the fact that the dimensional object a_B/" 2£ is invariant under the renormal- 
ization group transformations. [— ea s + /3(a s )] is the L>-dimensional /3-function, Z^J is the 
coefficient of the first e-pole in Z as defined below. 

Renormalization constants within the MS-scheme do not depend on dimensional param- 
eters (masses, momenta) || and have the following structure: 

= 1 + E^P» (6) 

n=l 

Since Z as does not depend explicitly on fj,, the /3-function is the same in all MS-like schemes, 
i.e. within the class of renormalization schemes which differ by the shift of the parameter \i. 
That is why the /3-function is the same in the MS-scheme ||] and in the MS-scheme ||. 

In general, the most straightforward way to obtain renormalization constants is by multi- 
plicative renormalization of the relevant Green functions 

rRenormalized(«s) = Z(a s ) Tsare^^), (7) 

This direct approach was used for the independent calculation of the 3-loop /3-function in Q. 
However the extension of this approach to 4-loops requires the calculation of 4-loop massless 
propagator type integrals which are still difficult to evaluate (at this moment). 
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Another approach to find the renormalization constants Zi is to obtain them as the sum 
of the counterterms of individual diagrams. This means that one applies the R-operation to 
each individual diagram that contributes to a propagator or vertex function. This allows a 
greater freedom in choosing the type of integrals that one needs to evaluate. This approach 
(based on massless propagator type integrals) was used for the first 3-loop calculation of the 
(5- function Q. 

For the calculation presented in this article we use massive integrals, instead of the mass- 
less integrals used in previous calculations. The calculation of renormalization constants 
within the MS-scheme can be reduced to the calculation of massive vacuum bubble integrals 
(i.e. massive integrals with no external momenta) using the general method of infrared rear- 
rangement [10]. This method uses the property that within dimensional regularization overall 
ultra-violet divergences are polynomial in external momenta and masses, also for individual 
diagrams. For the renormalization constants of the ghost-ghost-gluon vertex, ghost propa- 
gator and gluon propagator this means, that we can safely apply Taylor expansions in the 
external ghost- and gluon momenta if we introduce a non-zero auxiliary mass M for all inter- 
nal propagators (also for the gluons). It is understood that this auxiliary mass serves only as 
an infrared cutoff parameter that is nullified after renormalization of the individual diagrams. 
For simplicity, we introduce the mass only in the denominators of the propagators, not in the 
numerators. The difference between the overall divergences of the diagrams with and without 
the mass is a term that is polynomial in this mass and vanishes when M is nullified. 

The procedure of renormalization with an auxiliary mass works well for individual dia- 
grams. However the introduction of the mass M in the gluon propagators spoils multiplicative 
renormalizability of the (massive) Green functions. 

Here we use an intermediate approach to renormalization in order to get the 4-loop coun- 
terterms for the sum of the diagrams. We compute poles in e of the corresponding 4-loop 
massive diagrams. But we do not renormalize each diagram separately. The subtraction of 
sub divergences is done for the whole sum of the 4-loop diagrams. This is done by means of 
adding to the sum of the 4-loop diagrams the sum of the necessary bare diagrams of 1,2 and 
3-loops with all vertices replaced by effective vertices and all propagators replaced by effective 
propagators. The effective vertices contain the necessary vertex renormalization constants up 
to the appropriate order in a s and similarly, the effective propagators contain the necessary 
propagator counterterms. The various vertex and propagator renormalization constants that 
are needed in the effective vertices and propagators are already known from the lower order 
massless calculations (we emphasize that they are mass independent) except for the overall uv 
divergence proportional to M 2 of the gluon propagator, which is only needed up to 3-loops. 

Special routines for the symbolic manipulations program FORM [jll|] were constructed 
to efficiently evaluate the 4-loop massive bubble integrals up to pole parts in e and corre- 
spondingly of the 3-loop massive bubbles up to finite parts. For the 4-loop integrals, we 
only needed to deal with two master bubble topologies, see fig. 1. The various vertex and 
propagator diagrams were generated by means of the diagram generator QGRAF |0|. For 
the present calculation we evaluated of the order of 50.000 4-loop diagrams. 
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Figure 1. The basic (master) vacuum bubble topologies up to 4 loops. It 
is understood that all lines have the same non-zero mass, M. There are two 
basic 4- loop topologies (the last one shown is non-planar). 



We obtained in this way the following result for the 4-loop beta function in the MS-scheme 
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Here C is the Riemann zeta-function (C 3 = 1.202056903 • • •)• [T a T a } i:j = C F 5 i:j and f acd f bcd = 
C A 5 ab are the Casimir operators of the fundamental and the adjoint representation of the Lie 
algebra. ti(T a T b ) = T F 5 ab is the trace normalization of the fundamental representation. N A 
is the number of generators of the group (i.e. the number of gluons) and nf is the number 
of quark flavours. We expressed the higher order group invariants in terms of contractions 
between the following fully symmetrical tensors: 



d abcd 



1. 



:Tr 



rj^a,rj-ibrj-iCrj-id rj-iarj-ibrj-idrj-iC _|_ rj^CL rj-iC rj~ib rj-id 



rparj-iCrj-idrpb _|_ rparpdrpbrpC _|_ rj^arpdrpCrpb 



(9) 
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d abcd = i Tr 
A 6 



^ja^jb ^jc^jd _|_ ^ja^jb^jd^jc _|_ ^ja^jc^jb^jd 



+c a c c c d c b + c a c d c b c c + c a c d c c c 6 



(10) 



where the matrices [C a ]& c = —if abc are the generators in the adjoint representation. The 
result of eq. (||) is valid for an arbitrary semi-simple compact Lie group. The result for QED 
(i.e. the group U(l)) is included in eq.(|) by substituting Ca = 0, d a A bcd = 0, Cp = 1, Tp = 1, 
d abcd = 1; tv a = 1. This result for QED agrees with the literature |l3| . A second independent 
check of eq. (||) is provided by the calculation [|14|] where the large-nj terms for the QCD 
beta- function were calculated in all orders of the coupling constant. Our terms agree with 

0- 

The result of eq.(||) is obtained in an arbitrary covariant gauge for the gluon field. This 
means that we keep the gauge parameter £ that appears in the gluon propagator i [— g^ u + — 
Oq^Q" ' I '(l 2 + ^ e )]/( ( ? 2 + * e ) as a f ree parameter in the calculations. The explicit cancellation of 
the gauge dependence in the /3-function gives an important check of the results. The results 
for individual diagrams that contribute to the (3 function also contain (apart from the constant 
£3) the constants £4, (5 and several other constants specific for massive vacuum integrals. The 
cancellation of these constants at various stages in the calculation provides additional checks 
of the result. 

For the standard normalization of the SU(iV) generators we find the following expressions 
for the colour factors 

1 ^ Ar ^ N 2 -l d a A bcd d a A bcd N 2 (N 2 + 36) 
T F = -, C A = N, C F = ^ ]r , — = , 

df cd d a A bcd _ N(N 2 + 6) df cd df cd _ N 4 - QN 2 + 18 
Na~ ~ 48 ' Na~ ~ 96N 2 

Substitution of these colour factors for N = 3 into eq.(||) yields the following result for QCD 
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f3 2 ~ 1428.50 - 279.611ra/ + 6.01852nf 

[3 3 » 29243.0 -6946.30^+405.089^ + 1.49931^ (12) 

We note that (3% is positive for all positive values of nj. 

It is interesting to compare our result with a recent prediction |l5[] for the 4-loop coeffi- 
cient of the QCD /3-function in the MS-scheme using Pade Approximants. For nj = 3 this 
prediction is within a factor 2 of the exact result (for nj = 5 this factor is about 9). 

Considering that one might want to use the results of equation || for different groups or 
representations we will also express its constants in a different way In general one can 
write for a representation R of a simple Lie group 

d abcd = h{R)d abcd + hA31^ ab 5 cd + 5 ac 5 bd + 5 ad 5 bc ) 

3 

in which the tensor d is now traceless. The only exceptions to this are the spinor representa- 
tions of SO(8) for which there are two fully symmetric traceless tensors with 4 indices. We 
will not consider this special case here. The normalization of the tensor d is fixed by the 
definition of 14(F). Contraction with S ab S cd gives 

' v ' y N R &h{R) 

= T R (C R -\c A ) (13) 
6 

with I 2 {R) = T R , I 2 {A) =T A = C A and N A T R = C R N R . We use here that N R and N A are 
the dimensions of the representation R and the adjoint representation respectively. For the 
adjoint representation we have that J 2 ,2(A) = 5C A /6. 
After this the following identities hold: 

jabcdjabcd jabcd jabcd o 

jabcd mbcd jabcdjabcd q 

d -^- = h{R)iM)-^— + jt—MR)J2M) (14) 

M A N A N A + 2 

jabcd ^abcd jabcdjabcd q 

These tfs have some nice properties. They are zero for all exceptional groups and for SU(3). 
In addition they are representation independent. Hence we need to give them only for the 
classical groups: 

d a b cd d abcd {su{N)) = N A (N A -3)(N A -8) 
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d abcd d abcd {SO{N)) 
d abcd d abcd {SP(N)) 



N A (N A - 1)(N A - 3) 

12{N A + 2) 
N A (N A - 1)(N A - 3) 



(15) 



192{N A + 2) 

For the overall normalization we define the constant 6 with the relation 

Tr \c a C p ] = bg 5 a/S 



(16) 



in which g is the dual Coxeter number. The factor b is related to the normalization factor a 
in the article by Cvitanovic [17]. The relation is b = 2a for the groups SU(iV) and SP(iV) and 
b = a for SO(iV). The canonical choice of b is 1 for all groups. One should however be aware 
of the fact that sometimes different choices are used, especially for the exceptional groups. 
Some values for the fundamental and adjoint representations are: 
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For all groups we have Cf = TpN A /Np. Ia{F) = b 2 for the classical groups and zero for 
the exceptional groups. For other representations one would have to obtain values for the 
quantities Tr, Nr and Ia(R). By comparing the above values with the equations |] and || 
one may observe that the choice of a different value for b corresponds to a redefinition of the 
coupling constant. 
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